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Earthquakes (EQs) are large-scale fracture phenomena in the Earth's heterogeneous crust. Frac- 
ture induced physical fields allow a real-time monitoring of damage evolution in materials during 
mechanical loading. Electromagnetic (EM) emissions in a wide frequency spectrum ranging from 
kHz to MHz are produced by opening cracks, which can be considered as the so-called precur- 
sors of general fracture. We emphasize that the MHz radiation appears earlier than the kHz in 
both laboratory and geophysical scale. An important challenge in this field of research is to dis- 
tinguish characteristic epochs in the evolution of precursory EM activity and identify them with 
the equivalent last stages in the EQ preparation process. Recently, we proposed the following two 
epochs/stages model: (i) The second epoch, which includes the finally emerged strong impulsive 
kHz EM emission is due to the fracture of the high strength large asperities that are distributed 
along the activated fault sustaining the system, (ii) The first epoch, which includes the initially 
emerged MHz EM radiation is thought to be due to the fracture of a highly heterogeneous system 
that surrounds the family of asperities. A catastrophic EQ of magnitude Mw = 6.3 occurred on 
06/04/2009 in central Italy. The majority of the damage occurred in the city of L'Aquila. Clear 
kHz - MHz EM anomalies have been detected prior to the L Aquila EQ. Herein, we investigate the 
seismogenic origin of the detected MHz anomaly. The analysis in terms of intermittent dynamics of 
critical fluctuations reveals that the candidate EM precursor: (i) can be described in analogy with 
a thermal continuous phase transition; (ii) has anti-persistent behaviour. These features suggest 
that the emerged candidate precursor could be triggered by microfractures in the highly disordered 
system that surrounded the backbone of asperities of the activated fault. We introduce a criterion 
for an underlying strong critical behavior. In this field of research the reproducibility of results is 
desirable: the results should be verify by a number of precursory MHz EM emissions. We refer to 
precursory MHz EM activities associated with nine significant EQs that occurred in Greece. We 
conclude that all the MHz EM precursors under study also can be described in analogy with a 
continuous second-order phase transition showing antipersistent behaviour. 

PACS numbers: 



I. INTRODUCTION 

EQs are large-scale fracture phenomena in the Earth's 
heterogeneous crust. A vital problem in material sci- 
ence and geophysics is the identification of precursors of 
macroscopic defects or shocks. In physics, the degree to 
which we can predict a phenomenon is often measured by 
how well we understand it. However, despite the large 
amount of experimental data and the considerable ef- 
fort that has been undertaken by the material scientists, 
many questions about fracture processes remain stand- 
ing. Especially, many aspects of EQ generation still es- 
cape our full understudying. This fact is reflected in the 
disappointing progress on short-term EQ prediction [1]. 

Fracture induced physical fields allow a real-time mon- 
itoring of damage evolution in materials during mechan- 
ical loading. EM emissions in a wide frequency spec- 
trum ranging from kHz to MHz are produced by opening 
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cracks, which can be considered as the so-called precur- 
sors of general fracture. These precursors are detectable 
both at a laboratory and a geophysical scale [1-12]. Our 
main observational tool is the monitoring of the fracture 
which occurs in the focal area before the final break-up 
by recording their kHz-MHz EM emissions. 

Since 1994, a station was installed at a mountainous 
site of Zante island (37.76<W- 20.76°.E) in the Ionian Sea 
(western Greece) with the following configuration: (i) six 
loop antennas detecting the three components (EW, NS, 
and vertical) of the variations of the magnetic field at 3 
kHz and 10 kHz respectively; (ii) two vertical A/2 electric 
dipoles detecting the electric field variations at 41 and 54 
MHz respectively and (iii) two Short Thin Wire Anten- 
nas (STWA) of 100 m length each, lying on the Earth's 
surface, detecting ultra low frequency (ULF) (< 1 Hz) 
anomalies, at EW and NS direction respectively. The 3 
kHz, 10 kHz, 41 MHz, and 54 MHz were selected in order 
to minimize the effects of the sources of man-made noise 
in the mountain area of the Zante island. All the EM 
time-series were sampled at 1 Hz. Such an experimental 
setup helps to specify not only whether or not a single 
EM anomaly is preseismic in itself, but also whether a 
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combination of such disturbances at different frequencies 
can be characterized as pre-seismic. 

Clear EM anomalies have been detected over periods 
ranging from a few days to a few hours prior to recent 
destructive EQs in Greece, and it is by now recognized 
that the pre-fracture kHz-MHz EM time series contain 
valuable information about the fracture EQ preparation 
process [2-11]. Being non-destructive, monitoring tech- 
niques based on these fracture-induced fields are a basis 
for a fundamental understanding of fracture mechanism 
and for developing consecutive models of rock/focal area 
behaviour [12]. An improved understanding of the kHz- 
MHz EM precursors, especially its physical basis, has 
direct implications for the study of EQ generation pro- 
cesses and EQ prediction research. There is a growing 
number of indications on possible coupling between iono- 
spheric anomalies and EQs [3]. The appearance of the 
ULF anomaly verifies that a physically powerful fracto- 
electrification, and thus the fracture process, has been 
extended up to the surface layer of the crust. 

We emphasize that the MHz radiation appears ear- 
lier than the kHz activity [1, 10]. This situation is in 
harmony with laboratory results [1] . An important chal- 
lenge in this field of research is to distinguish character- 
istic epochs in the evolution of precursory EM activity 
and identify them with the equivalent last stages in the 
EQ preparation process. Recently, we have proposed the 
following two epochs/stages model of EQ generation [5]: 

(i) The abrupt emergence of strong impulsive kHz EM 
emission in the tail of the precursory EM activity is trig- 
gered by the fracture of the backbone of high strength 
and large asperities that are distributed along the acti- 
vated fault sustaining the system. The analysis of the 
detected precursory kHz EM activities indicate an un- 
derlying non-equilibrium process without any footprint 
of an equilibrium second-order phase transition. Notice, 
these precursory time series show persistent behavior, 
i.e., if the amplitude of EM fluctuations increases in a 
time interval, it is likely to continue increasing in the in- 
terval immediately following. The persistency implies a 
nonlinear positive feedback in the fracto-EM mechanism. 
The existence of a positive feedback justifies the absent 
of any signature of a second order transition. Based on 
well documented scaling properties of fault surface topol- 
ogy, we argue that the candidate kHz EM precursors 
might be originated during the slipping of two rough and 
rigid fractional-Brownian-motion-type profiles one over 
the other, with a roughness which is consistent with field 
and laboratory studies [7]. The results also suggest that 
the activation of a single EQ (fault) is a reduced self- 
affinc image of the whole regional seismicity and a magni- 
fied self-affine image of the laboratory seismicity [7] . The 
last suggestion is also supported by arguments based on 
non-extensive (Tsallis) statistical mechanics [8,9]. 

(ii) The initial MHz anomaly is thought to be due to 
the fracture of a highly heterogeneous system that sur- 
rounds the family of asperities [7]. The associated ar- 
guments have been supported by means of criticality, in 



particular in terms of the recently introduced method of 
critical fluctuations and Levy statistics [11]. More pre- 
cisely, the analysis indicates that the MHz EM activity 
is well described in analogy with a thermal continuous 
phase transition [5] . A fractal spectral analysis by means 
of Hurst exponent shows that this activity has antiper- 
sistent properties, i.e., if the EM fluctuations increase in 
period, they are likely to continue decreasing in the inter- 
val immediately following and vice versa [5, 10]. An an- 
tipersistent behavior implies a set of fluctuations tending 
to induce a stability within the system, i.e., the existence 
of a nonlinear negative feedback in the underlying fracto- 
EM mechanism which "kicks" the opening cracks away 
from extremes. The observed antipersistency, which is 
rooted in the heterogeneity of the system [5, 10], enables 
the fracture in highly heterogeneous systems to be de- 
scribed via an analogy with thermal continuous phase 
transitions. We have focused especially on the naturally 
arising question: what is the physical mechanism that 
organizes the heterogeneous system in its critical state? 
Combining ideas of Levy statistics, Tsallis statistics, and 
criticality on one hand and features hidden in the precur- 
sory MHz time series on the other hand, we argue that 
a Levy walk type mechanism can organize the heteroge- 
neous system to criticality [11]. Moreover, based on a 
numerically produced truncated Levy walk, we propose 
[11] a way to estimate in the stage of critical fluctuations: 
(i) the associated Levy index-a, which describes quantita- 
tively the underlying Levy dynamics, and (ii) the range 
of values where the nonextesitive Tsallis index q is re- 
stricted [11]. 

A catastrophic EQ of magnitude Mw = 6.3 occurred 
on April 6, 2009 (Olh 32m 41s UTC) in central Italy. The 
majority of the damage occurred in the city of L Aquila. 
Notice, this event was very shallow, namely, its depth was 
a few kilometres. Clear EM anomalies from ULF, kHz 
up to MHz have been detected prior to the L' Aquila EQ. 
Herein, we focus on the MHz EM anomaly that emerged 
on March 26 1999. The main goal of this contribution 
is to investigate the possible seismogenic origin of the 
detected anomaly in terms of the above mentioned two 
stages / epochs model. More precisely, we attempt to 
verify that the detected MHz EM activity: (i) is well 
described in analogy with a thermal continuous phase 
transition, and (ii) is characterized by Hurst exponents 
that indicate the existence of an underlying antipersistent 
fracture mechanism. The possible appearance of these 
two features taken together support the hypothesis that 
the recorded precursor coud be rooted in the fracture of 
the heterogeneous medium that surrounded the family 
of strong asperities that were distributed along the acti- 
vated fault of the L'Aquila EQ. The analysis reveals that 
the aforementioned two expected features are hidden in 
the precursor under study. 

The physics of EQs has been demonstrated to be a 
very complicated matter, and the road seems to be long 
and hard. A difficulty stems from the lack of large and 
reliable database: most studies in this area have been 
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limited by a lack of enough experimental results to con- 
duct a statistically significant analysis of the phenomena 
and obtain firm results. In this field of research the repro- 
ducibility of results is really desirable: the results should 
be verified by a numerous of precursory EM emissions. 
However, the collection of such a volume of appropriate 
EM data for statistical purposes requires some decades 
of years at least. Indeed, due to their absorption the EM 
precursors are associated with surface EQs of magnitude 
approximately 6 or greater that occur on land or closed to 
coastline . In the present study we investigate candidate 
precursory MHz EM activities associated with other 9 
EQs of magnitude approximately 6 or greater that occur 
on land or closed to coastline in Greece. The results are 
in surprising agreement with those found in the study of 
the precursor associated with the L'Aquila EQ. Indeed, 
all the MHz EM precursors can be described in terms 
of a continuous second order phase transition and show 
antipersistent behavior. 

We emphasize that kHz EM anomalies have also 
recorded prior to the L'Aquila EQ on April 4, 2009. The 
performed multidisciplinary analysis supports the con- 
sideration that this precursory EM activity could be trig- 
gered by the fracture of a backbone of high strength and 
large asperities that are distributed along an activated 
fault. Finally, the appearance of a clear ULF anomaly 
during the last days prior to the L'Aquila EQ strongly 
verifies that a physically powerful fracto-electrification, 
and thus a fracture process, has been extended up to the 
surface layer of the crust. We recall that the L'Aquila 
EQ was very shallow. The successive launch of ULF, 
kHz and MHz EM anomalies during a few days prior to 
the L'Aquila EQ helps to specify not only whether or 
not a single EM anomaly is preseismic in itself, but also 
whether a combination of such disturbances at different 
frequencies could be characterized as pre-seismic. A rel- 
evant manuscript is under preparation. 

The organization of this work is the following : In sec- 
tion 2 the method of critical fluctuations (MCF) [13, 14] 
is presented. In section 3 we apply the MCF to the de- 
tected MHz EM animaly prior to the L'Aquila EQ. In 
section 4 the reproducibility of results is investigated. In 
section 5 the connection between the antipersistency and 
criticality is demonstrated. Finally, in section 6 we dis- 
cuss our conclusions. 



II. THE METHOD OF CRITICAL 
FLUCTUATIONS 

The MCF, which constitutes a statistical method of 
analysis for the critical fluctuations in systems that un- 
dergo a continuous phase transition at equilibrium, has 
been introduced, recently [13, 14]. The authors have 
shown that the fluctuations <j) of the order parameter 
that corresponds to successive configurations of critical 
systems at equilibrium obey a dynamical law of intermit- 
tency, which can be described in terms of a 1-D nonlinear 



map. The invariant density p(4>) for such a map is charac- 
terized by a plateau which decays in a super-exponential 
way (see Fig. 1 in [13]). Importantly, the exact dynamics 
at the critical point can be determined analytically for a 
large class of critical systems introducing the so-called 
critical map [13]. For small values of (f>, this critical map 
can be approximated as 

<j>n+l = <Pn + U(pn + e n . (1) 

The shift parameter e„ introduces a non-universal 
stochastic noise: each physical system has its character- 
istic "noise" , which is expressed through the shift pa- 
rameter e n . Notice, for thermal systems the exponent z 
is connected with the isothermal critical exponent 5 as 
z = 5 + 1 . We emphasize that the plateau region of the 
invariant density p(<j>) corresponds to the laminar region 
of the critical map, where fully correlated dynamics takes 
place. The laminar region ends when the second nonlin- 
ear term in Eq. (1) becomes relevant. However, due to 
the fact that the dynamical law (1) changes continuously 
with <f>, the end of the laminar region cannot be easily 
defined based on a strictly quantitative criterion. Thus, 
the end of the laminar region should be generally treated 
as a varying parameter. 

Based on the aforementioned description of the criti- 
cal fluctuations, the MCF develops an algorithm permit- 
ting the extraction of the critical fluctuations, if any, in 
a measured time series. The crucial observation in this 
approach is the fact that the distribution P(l) of the lam- 
inar lengths I of the above mentioned intermittent map 

(1) in the limit e n — > is given by the power law [14] 

p(0~r«, (2) 

where the exponent pi is connected with the exponent 
z via pi = 

Therefore the exponent pi is connected with the 
isothermal exponent 5 as follows: 

Pi = 1 + ] (3) 

Inversely, the existence of a power law such as in Eq. 

(2) , accompanied by a plateau form of the correspond- 
ing density p{4>), is a signature of underlying correlated 
dynamics similar to the critical behavior [13, 14]. 

We give emphasis to the fact that it is possible in the 
frame of universality, which is characteristic for continu- 
ous second-order transitions, to give meaning to the ex- 
ponent pi beyond the thermal phase transitions [14]. 

The MCF is directly applied to time series or to seg- 
ments of time series which appear a cumulative type sta- 
tionary behavior. The main aim of the MCF is to esti- 
mate the exponent p/. The distribution of the laminar 
lengths of fluctuations included in a stationary window 
is fitted by the function: 
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P(l) ~ r^e-vzi ( 4 ) 

We focus on the exponents P2 and p$. If the exponent 
Pz is zero, then, the exponent pi is equal to the expo- 
nent pi. Practically, as the exponent ps approaches to 
zero, then, the exponent P2 approaches to pi, while, the 
laminar lengths tend to follow a power-law type distribu- 
tion. We note that Eq. (3) suggests that the exponent 
Pi ( OT P2) should be greater than 1. In conclusion, the 
critical profile of the temporal fluctuations is restored in 
the restrictions: p 2 > 1 and p 3 w 0. 

We stress that when the exponent p 2 is smaller than 
one, then, independently of the p3-value, the system is 
not at a critical state. Generally, the exponents p 2 ,P3 
have a competitive character, namely, when the expo- 
nent P2 decreases the associated exponent p% increases 
(they are mirror image each to other). In this way, we 
can identify the deviation from the critical state. To be 
more concrete, as the exponent p 2 (P2 < 1) is close to 1 
and simultaneously the exponent p 3 is close to zero, then, 
the system is at a sub-critical state. As the system re- 
moves from the critical state, the exponent p 2 decreases 
while simultaneously the exponent p^ increases reinforc- 
ing, in this way, the exponential character of the laminar 
lengths distribution. A situation like this indicates that 
the dynamics of the system becomes rather chaotic than 
deterministic [15, 16]. 

In summary, the research of criticality in natural sys- 
tems can be quantitatively accomplished by estimating 
the values of only two parameters, namely the exponents 
p 2 and p 3 . 

Up to now, the MCF has been applied on numerical 
experiments of thermal systems (Ising models) [14], on 
electromagnetic pre-seismic signals [5], and on electro- 
cardiac signals from biological tissues [17]. 

Comment: The end of laminar regions is a variable pa- 
rameter inside the decay region. In a recent work [18], it 
has been presented a prototype intermittent map of the 
form (1) which was in agreement with the provided prop- 
erties. In the frame of this prototype, the exponent p 2 
had almost constant value for all the end points of lami- 
nar regions. This means that every trajectory carry the 
correct information of the nonlinear term. This property 
could be characterized as a completeness behavior which 
notices a stable character of criticality. This behavior is 
rarely observed in real time series and so, we consider as 
p 2 exponent the maximum value because for this ^2-value 
the exponent p 3 has its minimum value and so the p 2 is 
closest to the exponent of criticality pi . 

Based on the above mentioned considerations, the case 
where the p 2 -values are greater than 1 and the corre- 
sponding p3-values close to zero for all end points of lam- 
inar regions indicates that all the trajectories are critical 
and carry almost the same information about the dy- 
namical term. We can say that in this case the system is 
characterized by a strong criticality. 



III. APPLICATION OF THE MCF METHOD TO 
THE DETECTED PRE-SEISMIC MHZ EM 
RADIATION PRIOR TO THE L'AQUILA 
EARTHQUAKE 

As it was said, a catastrophic EQ of magnitude Mw 
= 6.3 occurred on April 6, 2009 (Olh 32m 41s UTC) in 
central Italy. The majority of the damage occurred in 
the city of L'Aquila. Notice, this event was very shal- 
low, namely, its depth was a few kilometres. Clear EM 
anomalies from ULF, kHz up to MHz have been detected 
prior to the L' Aquila EQ. Herein, we focus on the 41- 
MHz EM anomaly that emerged prior to the L'Aquila 
EQ, i.e., on March 26 2009. The main goal of this con- 
tribution is to investigate the possible seismogenic origin 
of the detected anomaly in terms of the above mentioned 
two stages / epochs model [5]. The data were sampled 
at 1 Hz. 

It is significant to isolate the window of the time series 
carrying the footprints of the critical state. Our analysis 
verified the existense of a critical window (CW). Indeed, 
we have identified a characteristic time interval including 
15000 points (Fig. la) in which a cumulative stationary 
condition occurs. As mentioned above, the existense of 
cumulative stationary behaviour in the candidate CW is 
a necessary condition for applying the MCF method to a 
nonequilibrium process. The stationary behaviour of the 
CW has been checked by estimating the mean value and 
the standard deviation for various time intervals in the 
CW, all having a common origin. The results are shown 
in Fig. (lb). The corresponding distribution of the am- 
plitude P(ip) of the emerged EM pulses in the candidate 
CW is shown in Fig. 2. It is characteristic the appearence 
of the plateau region in the top of distribution, as it is 
provided for the invariant density of critical map [13]. 
The next step is to produce the distribution P{1) where 
the length 1 is the laminar interval, i.e., the stay time 
within the laminar region. More precisely, the laminar 
lengths are described by the lengths of the sub-sequences 
in the EM time series that are resulted from successive 
^-values obeying the condition ip < ip < ipi, where ip Q 
is the fixed point and the end of the laminar region ipi 
is a variable parameter. The more abrupt decay side of 
the distribution P(ip) is usally the fixed point i\> [15] but 
in the case of a symmetric distribution, like here, we se- 
lect as fixed point the beggining of the distribution (here 
ipo = 650). Using the fitting function (4) we estimate 
the exponents P2,P3 for different -0;- values. In Fig. 3a 
the exponents P2,P3 vs the end point ipi are shown. In 
order to avoid end effects we ignore the exponents values 
coming from the end points of the distribution P(ip)- As 
it results from Fig. 3, this window can be characterized 
as a strong CW. Indeed, the majority of the j^-values are 
restricted in a zone from 1.20 to 1.43, while the exponent 
Ps has values close to zero. Notice, that the quality of 
fitting in fig.3b was excellent (R 2 > 0.99). 

In fig. 4, for comparison reasons, we present a typical 
distribution P(Z) for a time interval (28/3/2009) regarded 
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FIG. 1: The case of the L'Aquila EQ. (a) The MHz EM pre- 
cursory fluctuations window . This window has a duration 
of 15000 points (lp/sec) and appeared almost 3.5 days before 
the EQ occurence, (b) The temporal evolution of the cumula- 
tive mean value and standard deviation for increasing number 
of points starting from an initial set of 6000 points. 




FIG. 2: The distribution P(tp) for the values of timeseries 
1(a) is shown . 



FIG. 3: (a) Case of the L'Aquila EQ. The exponents P2,P3 
vs ipi indicate an underlying strong critical behavior. The 
corresponding critical window has a length of 15000 points 
(lp/sec) and is appeared almost 3.5 days before the EQ oc- 
curence. The p2-values are restricted in the zone (1.20-1.43). 
(b) The distribution of the laminar lengths, as it is produced 
by the analysis in terms of MCF, for end point ipi = 770. The 
fitting function gives p2 = 1.43 and P3 = 0.0008. For concrete 
end point the distribution P(l) is more close to a power-law 
(p3 minimum and p2 maximum ). The solid line represents 
the fit of data by the fitting function (Eq. 4). The quality of 
fitting is very nice (R 2 = 0.996). 




as noise which does not exhibit criticality. 



In summary, based on the above mentioned results in 
terms of MCF, we argue that the candidate critical win- 
dow under study in the preseismic time series really is 
a critical window reflecting an underlying fracture pro- 
cess which is characterized by "strong criticality". The 
underlying dynamics posses intermittent characteristics, 
shown to describe the temporal fluctuations of the order 
parameter of a macroscopic system at its critical point. 



FIG. 4: The distribution of the laminar lengths, as it is pro- 
duced by the analysis in terms of MCF, from a segment of 
preseismic EM signal of L'Aquila (28/3/2009) which does not 
exhibit criticality. The fitting function gives p2 = 0.48 and 
P3 = 0.22. 
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IV. THE REPRODUCIBILITY OF RESULTS 

In this field of research the reproducibility of results 
is desirable: the results should be verified by a num- 
ber of precursory EM emissions. Thus, a question ef- 
fortlessly arises is whether the appearance of windows 
showing strong criticality, as it has been described previ- 
ously is systematically observed in candidate precursory 
MHz EM activities associated with significant EQs that 
occurred in Greece during the last years. In the present 
section we refer to MHz EM activities recorded prior to 
nine significant EQs that occurred in land or near the 
coast line of Greece. We note that it is possible to have 
the existence of more than one critical window in the 
preseismic EM fluctuations. In the following, for statis- 
tical reasons, we present the critical window which has 
the greater duration. 

(a) The Athens EQ. On 7 September 1999 (11:56:49 
UT) the Athens EQ (38.12°N, 23.60°E) occurred with 
magnitude M — 5.9 [5, 7,10]. In figs. 5a, 6a the associated 
CW has a duration 5.5h and it is appeared almost 4 days 
before the EQ event. 

(b) The Kozani-GrevenaEQ. On 13 May 1995 (8:47:13 
UT) the Kozani-Grevena EQ (40.17°N, 21.68°E) oc- 
curred with magnitude M = 6.6. In figs. 5b, 6b the 
associated CW has a duration 8h and it is appeared al- 
most 2 days before the EQ event. 

(c) The Zante EQ. On 11 April 2006 (17:29:28 UT) the 
Zante EQ (37.68°N, 20.91°E) occurred close to the Zante 
island with magnitude 5.9. In figs. 5c, 6c the associated 
CW has a duration 5h and it is appeared almost 19h 
before the EQ event. 

(d) The Mcthoni EQ. On 14 February 2008 (10:09:23 
UT) the Methoni EQ (36.55°N, 21.77°E) with a magni- 
tude M = 6.9 occurred. In figs. 5d,6d the associated 
CW has a duration 1.5d and it is appeared almost 6 days 
before the EQ event. 

(e) The Andravida EQ. On 8 June 2008 (12:25:28 UT) 
the Andravida EQ (37.95°N, 21.53°E) with a magnitude 
M = 6.4 occurred. In figs. 5e,6e the associated CW has 
a duration 5h and it is appeared almost 2 days before the 
EQ event. 

(f) The Lamia EQ. On 13 December 2008 (11:27:19 
UT) the Lamia EQ (38.72°N, 22.5°E) with a magnitude 
M — 5.7 occurred. In figs. 5f,6f the associated CW has 
a duration 6h and it is appeared almost 20h before the 
EQ event. 

(g) The Kefalonia EQ. On 25 March 2007 (13:57:58 
UT) the Kefalonia EQ (38.34°N, 20.42°E) with a magni- 
tude M — 6 occurred. In figs. 5g,6g the associated CW 
has a duration 5h and it is appeared almost 3 days before 
the EQ event. 

(h) The Agios Nikolaos EQ. On 13 Junuary 2009 
(06:12:43 UT) the Agios Nikolaos EQ (35.66°N, 26.39°E) 
with a magnitude M — 5.7 occurred. In figs. 5h,6h the 
associated CW has a duration 8h and it is appeared al- 
most 5 days before the EQ event. 

(i) The Strofades EQ. On 16 February 2009 (23:16:38 
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FIG. 5: The exponents P2,P3 vs ipi indicate an underlying 
strong critical behavior for the other 9 EQ events as these are 
described in text. The corresponding CW are the greater in 
duration for each EQ. 



UT) the Strofades EQ (37.13°N, 20.78°E) with a magni- 
tude M — 6 occurred. In figs. 5i,6i the associated CW 
has a duration 5.5h and it is appeared almost 10 days 
before the EQ event. 

We clarify that the precursors associated with the 
above mentioned (a)-(f) and (i) seismic events have been 
recorded by the sensors of the Zante station, while those 
associated with the (g) and (h) shocks have been de- 
tected by the sensors of the Neapolis (east Crete) sta- 
tion (35.26°N, 25.61°E) and Kefalonia station ((38.18°N, 
20.59°E), respectivelly. 

Figs. 5 and 6 importantly show that in all the cases 
the EM fluctuations included within the critical windows 
are fitted by the function (4) with exponent P2 having 
values greater than 1 and exponent pz having values 
close to zero, while the quality of fittings were excelent 
(i? 2 > 0.99). Based on this result, we conclude that in 
the EM precursor associated with the L'Aquila EQ are 
hidden features of "strong criticality" similar to those 
hidden in EM precursors associated with other nine seis- 
mic events. Thus the reproducibility of results is satisfied 
in ten different MHz EM precursors. 



V. THE PERSISTENT BEHAVIOUR OF THE 
CANDIDATE ELECTROMAGNETIC 
PRECURSOR 

The Hurst exponent H characterizes the 
persistent/anti-persistcnt properties of the signal 
[19-20]. The range < H < 0.5 indicates an anti- 
persistency reflecting that increases in the value of a 
time series are likely to be followed by decreases and 
conversely. Physically, this implies a set of fluctuations 
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FIG. 6: The distribution of the laminar lengths, as it is pro- 
duced by the analysis in terms of MCF, from the other 9 
EQ events assosiated to the corresponding critical windows 
of Fig. 5. All the distributions are more close to power-law, 
for each case. 

tending to induce a stability to the system, namely a 
non-linear feedback mechanism that "kicks" the opening 
rate away from extremes. On the contrary, the range 
0.5 < H < 1 indicates persistent behavior. This means 
that increases in the value of a time series are likely to 
be followed by increase, namely the system has been 
starting to self-organize by a positive feedback process. 

In Rcf. [5] we have introduced the following connection 
between the exponents H and p 2 : 

H = 2- \p 2 . (5) 

Eq. (5) suggests that the allowed range of p 2 values 
for a second-order phase transition, i.e., p 2 > 1, leads 
to the condition H < 0.5 which indicates an underlying 
antipersistent mechanism. On the other hand, Eq. (5), 
as it has been found [5] is valid only if the exponent p 2 
lies in the interval 1 < p 2 < 1.5. The present analysis 
systematically leads to exponents p 2 which arc in agree- 
ment with the above limits of p 2 for all the 10 EQs under 
study. 

To summarize, in the EM precursor possibly associ- 
ated with the L'Aquila EQ are hidden features of strong 
criticality and antipersistency similar to those hidden in 
the EM precursors associated with other 9 seismic events. 
Notice, the interplay between the heterogeneities in the 
pre-focal area and the stress field could be responsible for 
the observed antipersistent pattern [5, 10 ]. 

VI. CONCLUSION-DISCUSSION 



terials during mechanical loading. These EM precur- 
sors are detectable both at a laboratory and a geological 
scale containing valuable information about the fracture 
/ EQ preparation process. MHz EM anomalies have been 
recorded prior to the catastrophic L'Aquila EQ a few 
days before the EQ occurrence. EM anomalies possibly 
associated with the preparation of EQs would be recog- 
nized as real precursors only when the physical mecha- 
nism of their origin would be clarified. Thus, an impor- 
tant challenge in this field of research is to connect the 
detected precursory EM activity with a last stage in the 
EQ preparation process. Previous studies have suggested 
that the systematically detected MHz EM activity prior 
to significant EQs could be originated during cracking 
in the highly heterogeneous material that surrounds the 
backbone of large and strong entities distributed along 
the activated fault sustaining the system. We find that 
the candidate EM precursor associated with the L'Aquila 
EQ is in harmony with this suggestion. Indeed, in anal- 
ogy to the study of critical phase transitions in statistical 
physics, it has been proposed that the fracture of het- 
erogeneous materials could be viewed as a critical phe- 
nomenon [21-22]. Based on a recently presented statis- 
tical method of analysis for the critical fluctuations in 
systems which undergo a continuous phase transition at 
equilibrium [13-14], we show [5] that the detected MHz 
anomaly prior to the L'Aquila EQ can be described in 
analogy with a thermal continuous phase transition by 
means of a recently introduced method of critical fluc- 
tuations. Especially, we introduce a criterion for an un- 
derlying strong critical behavior. According to this crite- 
rion the majority of trajectories in the properly defined 
laminar region, "carries" the information of an underly- 
ing criticality. The precursor under study follows this 
criterion. Importantly, the precursor shows antipersis- 
tent behaviour indicating an underlying non-linear feed- 
back system that "kicks" the opening rate of cracks away 
from extremes. The antipersistent behavior is in agree- 
ment with a system which undergo a continuous phase 
transition at equilibrium. In this field of research the 
reproducibility of results is desirable. We study precur- 
sory MHz EM activities associated with nine significant 
EQs that occurred in Greece during the last years. We 
show that all these nine MHz EM precursors behave as 
that detected prior to the L'Aquila EQ: they can be de- 
scribed in analogy with a continuous second-order phase 
transition showing strong criticality and antipersistent 
behaviour. The thesis underlying the present effort is 
that the emerged MHz EM anomaly prior the L'Aquila 
EQ could be emitted by the fracture of disordered ma- 
terial in the focal area of this EQ that surrounded the 
family of asperities of the activated fault. 



Fracture induced EM fields ranging from kHz to MHz 
allow a real-time monitoring of damage evolution in ma- 
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